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Summary 

Earth-orbiting spacecraft often contain solar arrays or 
antennas supported by a preloaded mast. Because of weight 
and cost considerations, the structures supporting the space- 
craft appendages are extremely light and flexible; therefore, it 
is vital to investigate the influence of all physical and struc- 
tural parameters that may influence the dynamic behavior of 
the overall structure. The study reported in this paper primar- 
ily focuses on the mast for the space station solar arrays, but 
the formulations and the techniques developed in this study 
apply to any large and flexible mast in zero gravity. Further- 
more, to determine the influence on the circular frequencies, 
the mass moment of inertia of the mast has been incorporated 
into the governing equation of motion for bending. A finite 
element technique (MSC/NASTRAN) has been used to verify 
the formulation derived in this paper. Results indicate that 
when the mast is relatively flexible and long, the mast mass 
moment of inertia influences the circular frequencies. 

Introduction 

With the evolution of man-tended or manned spacecraft 
such as Space Station Freedom comes the need for generating 
more electrical power than for an ordinary satellite. More 
electrical power simply translates to having a greater number 
of solar arrays with larger surface areas. In order to have a 
large surface area solar array and yet be able to package the 
folded solar array in a launch vehicle, the solar arrays should 
be designed with a relatively large aspect ratio. Because solar 
arrays for manned spacecraft require a large aspect ratio, their 
supporting structures (preloaded masts) will be as tall as and 
much more flexible than the mast for an ordinary satellite. 
Therefore, the assumptions made to formulate the governing 
equation of motion for a compact solar array will not be 
applicable for a solar array with a large aspect ratio. 
Although extensive studies have been performed on 
preloaded mast vibrations (refs. 1 and 2), the influence of the 
mast mass moment of inertia has not been considered. 

Our studies focused on developing the governing equation 
of motion for a preloaded beam, representing a typical mast, 
with a tip mass, representing the storage boxes of the solar 
cells (fig. 1.) The equation of motion for the beam was used 
to derive the eigen condition (characteristic equation) whose 
roots are the circular frequencies of the structure in question. 



Figure 1 . — Space station split blanket solar array. 


The accuracy of these circular frequencies is the key in deter- 
mining the forces and moments resulting from an outside 
excitation, such as space shuttle docking loads. The loads 
determined from the circular frequencies with high modal 
mass participation are often used to design the structure. 

Symbols 

C constant 
EA axial stiffness 
EI CS flexural stiffness 
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F force 

/ moment of inertia of mast cross section 

K preload parameter 

A '2 tip mass preload parameter 

L length of mast 

M moment 

Mj tip mass 

m mass of mast 

m' distributed mass 

P preload 

P[> component of preload in lateral direction 
Q shear load 
T function of time 
t time 

V lateral deformation as function of beam length and time 

v lateral deformation as function of beam length 

x distance 
P frequency parameter 
S axial deflection 
O slope of deflection 
X characteristic value 

p mass density of the mast 

(O circular frequency 
d partial derivative 

Subscript: 
n rtth root 

Superscripts: 

' differentiation with respect to length of mast 
. differentiation with respect to time 

Theory 

The solar arrays masts are essentially bars that are subject 
to bending vibrations. The source of the vibrations may be the 
interaction of the propulsion system, and/or the shuttle docking. 
Since a bar is an elastic body, its mass and stiffness are distrib- 
uted continuously along its length. For an elastic body the 
number of the generalized coordinates needed to describe the 
motion of the bar becomes infinite, and, as a result, the number 
of modes and mode shapes become infinite. 

To formulate the equation of motion for and analyze a bar 
under a bending vibration, the following assumptions were 
made: (1) the material is homogeneous and isotropic, and it 
follows Hooke’s law; (2) the continuum representing the mast 
is conservative, and the structural damping and viscous 


damping (frictional forces) are nearly zero; (3) the magnitude 
of displacement is small when compared to the mast length, 
so the sine and tangent of the angle arc equal to the angle of 
deflection; (4) the mast has a constant cross-sectional area 
and constant bending stiffness along its length; and (5) the 
supporting spacecraft is rigid and massive when compared to 
the solar array mast, and the mast is fixed at its base. It is 
important to note that the boundary conditions for the mast 
are derived based on the last assumption. 

The mechanism for preloading the mast (fig. 1) can be 
described as the tensioning of the solar array substrates by 
some force at their edges; a compressive force of equal mag- 
nitude will act on the mast. The substrates are tensioned to 
achieve flat surfaces and to increase their natural frequency to 
some nonzero value. 

In this study the assumption has been made that the 
masses of the solar cell storage boxes are concentrated at only 
one point, the tip of the mast. Such an assumption is appro- 
priate when analyzing the mast under a bending vibration 
since the storage boxes will experience the same direction 
and magnitude of motion during any bending mode. In con- 
trast, one tip mass is not a valid assumption when analyzing 
the mast under a torsional vibration, and an appropriate mass 
moment inertia of the solar cell storage boxes should be con- 
sidered. In summary, the continuum representing the mast is 
a preloaded mast that is fixed at one end and free at the other, 
where the tip mass is placed. 



Governing Equation of Motion 

Figure 2 shows a mast of length L, whose cross-sectional 
area is A and bending stiffness is EI CS . The preload has a mag- 
nitude P in the direction shown, and the tip mass has a mag- 
nitude Mj. One very important point to be made here is that 
because of the mechanism for preloading the mast, the direc- 
tion of the load always remains towards the root of the mast, 
regardless of the deflection and the mast mode shapes. This 
occurs because the solar cell substrates are essentially mem- 
branes and they react only to tensile loading. The magnitude 
of the preload remains constant, regardless of the magnitude 
of deflection, since the solar cell substrates contain negator 
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Figure 3. — Free-body diagram of a differential element of the 
preloaded beam. 


spring mechanisms at their edges that keep the magnitude of 

the preload constant. 

Figure 3 shows the free-body diagram of a differential 
element of the preloaded mast with all the forces and 
moments acting on its body. The sum of the loads can be 
expressed as follows: 


d 2 \ 


r , <? 4 V „<? 2 V , d 

c 'dx 4 dx 2 P dt 2 [dx 2 


= -pA 


d 2 \ 

dt 2 


(7) 


Since the displacement V is a function of length x and time t, 
the method of separation of variables is used: 

V = f(x,t) = \(x)T(t) (8) 

For a harmonic excitation 

7(0= sin cur (9) 

where a> is the circular frequency of motion. 

Substituting equation (9) into equation (8) and then differ- 
entiating with respect to x and t and substituting the results 
into equation (7) gives the following: 


'v + plco 2 '' 

v"-{ pA( ° 2 ) 

K Ef C s ) 

( EIcs J 


( 10 ) 


IF =0, which results in P = P 


ZF V =0, 


which results in 


d 2 \ 

-3x~- pA TF 


dQ 


0 ) 

( 2 ) 


Equation (10) is similar to equation (7), but it represents the 
motion at any given time, since v is a function only of x. 


XM o =0, 


which results in Q = 


<?M n dV 

+ P — 

dx dx 


-plQ 


(3) 


Using the elementary theory of strength of material gives 


and 


M = EI CS 


d 2 \ 

dx 2 



(4) 


(5) 


By differentiating equation (4) with respect to x once, 
differentiating equation (5) with respect to t twice, and then 
substituting the results into equation (3), we obtain 


Q = £/„^ + P^-p/4l 


f <? 2 v) 


<9x J 


dx 


dx 


dt 2 


( 6 ) 


Finally, differentiating equation (6) with respect to x and sub- 
stituting it into equation (2) gives the equation of motion for 
a preloaded mast with a tip mass: 


Discussion of Terms in Equation (10) 


When the preload and the mass moment of inertia are 
zero, equation (10) reduces to the equation of motion for a 
simple beam. In equation (10) the first term, v-9, represents 
the bending motion of a differential clement of length dx; the 


third term, 


f p Act ) 2 ^ 
V EI cs ) 


i V , represents the linear inertia of a differ- 


ential element of length dx that opposes the motion; and the 


second term. 


P + pIco‘ 
El 


.2 \ 


, is the contribution of the preload 


and the mass moment of inertia. The physical meaning of the 
second term can be described as follows: (1) When the mast 
is deformed laterally, the preload P tends to deform (bend) it 
even further, and the magnitude of this added displacement is 
proportional to the preload times the displacement lever arm; 
and (2) when the mast is deformed laterally, the angular 
acceleration tends to assist the bending of the mast in the 
same direction as the preload. It is interesting to note that the 
magnitude of this bending effect is influenced by the dis- 
placement v and by the circular frequency at which the sys- 
tem is vibrating. For high circular frequencies, this effect 
dominates the influence of preload. 
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General Solution 

For simplicity, let us define 


K 2 


P + p/ft > 2 


El 


cs 


and 


4 = pAco 
P EI„ 


(ii) 


( 12 ) 


Preload 



Figure 4. — Deflected shape of the preloaded beam centerline. 


where j3 is the frequency parameter (ref. 2). 

Substituting equations (11) and (12) into equation (10) 
results in 

v""+ K 2 \"-P 4 v = 0 (13) 


the free end is required. Figure 4 shows the mast as loaded; in 
this figure, Pd is a component of preload whose magnitude 
can be determined as follows: 


Equation (13) is a fourth order ordinary differential equation 
with the general solution 


tan 


L-8 


(19) 


v = Cj cosh AjX + C 2 sinh AjX 

+ C3 cos A 2 x + C 4 sin A 2 x (14) 


and 


P D = P sin 0 


( 20 ) 


where 


and 



Note that when the mast is deflected, its length is shortened 
by some magnitude 5. 

For a small deflection, 0 = tan 0 = sin 0; therefore 


P D = P sin 


0 = P0 = P 


V(L) 

L-8 


( 21 ) 


h=]i 2 


(16) 


The constants C\ to C4 can be determined from the bound- 
ary conditions and the initial condition. 

At the fixed end of the mast, the displacement and the 
slope of the displacement are zero; thus the first two bound- 
ary conditions are 


v(0) = 0 

(17) 

* =0 
dx *=« 

(18) 


Because the boundary condition at the free end, where the 
tip mass is placed, is more complex, a free-body diagram of 


M t V(L) 



t 

P D 

Figure 5. — Free-body diagram at free 
end of the preloaded beam. 
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where 5 is the axial deflection whose magnitude is influenced 
by the axial preload. Hs magnitude (ref. 3) is 

p i 

5 = — ( 22 ) 

EA 

Substituting equation (22) into (21) gives 


C/|^A 3 + A] A 2 jsinh AjL + ^A 2 + A 2 — A 2 j sin A 2 L 
+ /C 2 (cosh Aj L - cos A 2 L) 


+ C 7 


(A? + A,A 2 2 ) cosh AjL 


Pn = P 


yw 

PL 

E4 


L- 


(23) 


+ (_A 2 A?+A 3 2 ) 


A] « 

— - cos A?L 

a 2 


Figure 5 depicts the free-body diagram at the mast’s free 
end. To eliminate the time function, the terms from equa- 
tions (8) and (9) should be substituted for the terms in this 
free-body diagram. Summing the forces and moments at the 
mast free end results in 


A 


A 

+ K 2 | sinh A|L-— *-sin A 2 L 

*■2 j 


= 0 


(29) 


Q(L) + 


P y(L) 

L-* 

EA 


M t (i) 2 \(L) = 0 


(24) 

where 


and 


K,= 


M r 0) ~j~ 

EI„ 


(30) 


M(L) = EI CS 


dx 2 


= 0 

*=L 


(25) 


Note that since the direction of the preload always remains 
towards the mast fixed end, the moment caused by the pre- 
load about the mast fixed end is zero. 

If the mast is tall and flexible so that EA » PL, then 
equation (24) can be simplified to 


Q(L) + 


P y(L) 

L 


- M t co 2 \(L) - 0 


(26) 


Note that the following relations, which can easily be 
checked, simplify equation (29): AjA 2 = P 2m , A 2 +/C 2 = A 2 ; and 
A 2 — K 2 = A 2 . Thus equation (29) becomes 

C, [Aj A 2 sinh AjL — A 2 A 2 sin A 2 L + LT 2 (cosh A ] L — cos A 2 L)] 




Aj A 2 cosh AjL + A 3 cos A 2 L 


Therefore, equations (25) and (26) are the boundary condi- 
tions to be used at the mast free end. Using the first two 
boundary conditions (eqs. (17) and (18)) results in 


+AT?( sinh AiL — sinh A 2 L 

I a 2 


= 0 


C 3 =-C, 


(27) 


(31) 


and 

Q=-T- C 2 ( 28 ) 

Employing the third condition (eq. (26)) and equation (6) 
gives 


Employing the fourth condition (eq. (25)) results in 
Uj | A 2 cosh AjL + A 2 cos A 2 l| 

+ C 2 (A 2 sinh AjL-f AjA 2 sin A 2 l) = 0 


(32) 
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Table I. - Parameters Which Define the 
Idealized Solar Array Structure 
(ref. 4] 


Length, in 1146 

Width, in 380 

Beam El , psi 2.85x 1 0* 

Beam, psi ........ 6.9x 10* 

Beam weight, lb ........ 286.5 

Blanket weight, lb 766 

Tip weight, lb 67 

Tip inertia, lb-in. 2 8.06x1 0 5 


Solving equation (29) for C 2 and substituting that into equa- 
tion (34) yields 


v = Ci 


(cosh AjX -cos A 2 x) 


Equations (31) and (32) represent an eigenvalue problem. 
The determinant of C\ and C 2 must be zero, and this con- 
dition establishes the characteristics equation (having an 
infinite number of roots) that determines the characteristic 
number or frequency value. Setting the determinant of equa- 
tions (31) and (32) to zero gives 


A| A 2 sinh A | L - A^A 2 sin A 2 L + A" 2 (cosh A,L - cos 


Aj A 2 cosh A | L + A~] cos k 2 L + K 2 


sinh AtL sin 



2 p 6 -K 2 /} 4 sinh A,L sin A 2 L 


+/3 2 ^K a + 2f5 A j cosh Aj L cos A 2 C - + A 2 j 

x(A] cosh A ]L sin A 2 L - A 2 sinh AjL cos A 2 C) = 0 


x 


A ^ 

sinh A»x -sin A 9 x 

a 2 J 


( 35 ) 


(33) 

Equation (33) is similar to the eigen condition derived in 
reference 2 , with some minor differences in definition of 
parameters. Note that all parameters in equation (33) are func- 
tions of (o, and employing the constants from table I will result 
in a highly nonlinear equation with an infinite number of roots. 


Mode Shapes 

To determine the mode shapes, the constants C\ to C 4 
should be evaluated and substituted into equation (14). The 
boundary conditions of the mast can provide the necessary 
information to evaluate any set of three constants, but the last 
constant can be evaluated only with the initial condition. Since 
we are interested in the mode shapes at any instant of time and 
since the last constant can have any value, depending on the 
initial condition, an arbitrary value of 1.0 for the last constant 
is assumed in this study. 

Substituting equations (27) and (28) into equation (14) 
results in 


Equation (35) is the equation of motion for individual circular 
frequencies. When a circular frequency obtained from equation 
(33) is substituted into equation (35), the mode shape for that 
particular frequency can be determined. Note that, as stated 
previously, the value for C] is assumed to be 1 .0 prior to mode 
shapes determination. 

When the mast is excited, it will vibrate through all frequen- 
cies; therefore its motion, which is the sum of all mode shapes, 
can be expressed as 



(cosh A lw x - cos A 2ff x) 


A| fl A, n sinh A in L - A^A^ sin A -, n L + K 2 (cosh A ln L - cos A 2fJ L) 


l 3 A j n 

A. cosh A. L + A. cos A 2 L + K 2 sinh A . Z, sin A^, L 


V = C| (cosh 


A|X-cos A 2 x) + C 2 


sinh A|X 




(34) 


/l. 

x sinhA tn x- sin A 2fJ x 


(36) 
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where the subscript rt refers to the nth circular frequency. 

Again, it should be emphasized that equation (36) is the 
mode shape at any instant of time; to obtain the mode shape 
as function of the mast length and time, use equations (8) and 
(9), which yield 



(cosh X irj \ - cos A 2fI x) 



Figure 6. — A representation of the finite element model of the 
preloaded beam. Note: The mast and single beam element 
are physically at the same location, but In this figure they have 
been shown separately to emphasize the physics of loading. 
Node 50 is a shared node with a hinge. 


X\„X 2n S * n ^ ^ \n ^ j' I n^~2n S * n ^2n ^ X\ COS ^2/f^) 


X, n X\ n cosh X U L + X] n cos X 2n l + K 2 


sinh X in L - — — sin X 2n l 


x 


sinh X ln \ - 


X } 

~r ~ sin 
Xi„ 


sin 0) n t (37) 


where C n is obtained by using the initial condition. 


Finite Element Approach 

Since the normal mode solution of MSC/NASTRAN 
(Version 65c) does not have the capability of predicting the 
normal modes of a preloaded structure in one computer run, 
the general solution was divided into two separate runs. First 
the stiffness matrix of the preloaded mast was generated by 
using the geometric nonlinearity of solution 64 of MSC/ 
NASTRAN, and then the superelemcnt normal mode, 
solution 63, was used to determine the modes. 


constraining all the nodes in the z-direction. Node 1, which 
represents the fixed end of the mast, is constrained in all 
directions. 

To generate the preload on the mast, and yet force the 
direction of the preload to remain towards the fixed end of 
the mast, a single BEAM element with a very large stiffness 
was generated, running the entire length of mast. The single 
BEAM element shares node 50 with the mast. The large stiff- 
ness of the single BEAM element ensures that no mode cou- 
pling will occur between the mast and the single BEAM 
element. The node at the other end of the single BEAM ele- 
ment is node 100, which has been located 0.094 in. above 
node 1, in the positive x-direction. The reason for this offset 
is that when a preload (see table 1) is applied to the mast, 
because of the axial flexibility of the mast, the entire assem- 
bly will deform by 0.094 in. towards the negative x-direction. 
To enforce the motion of node 100 in the direction of the 
preload, a constraint in the v-direction was applied to this 
node. In addition, to eliminate the bending action between the 
single BEAM element and the mast, a hinge was placed at 
node 50. 

Finally, the database containing the preloaded stiffness 
matrix of the mast was used in solution 63 to obtain the circu- 
lar frequencies. The Generalized Dynamic Reduction tech- 
nique was used to extract the circular frequencies. 


Static Nonlinear and Normal 
Mode Solutions 

The finite element model of the mast, shown in figure 6, is 
approximately a 100 degree-of-freedom model, which con- 
tains BEAM elements. BEAM elements are .two- nod ed 
elements with the capability of deforming axially, laterally, 
and torsionally. For this study the torsional deformation capa- 
bility of the finite element model was eliminated by con- 
straining all the nodes from rotating about the x-direction. 
Also the model was forced to bend in only one direction by 


Results 

A Fortran program was developed to solve equation (33) 
for circular frequencies by using the Newton method. The 
circular frequencies from the Fortran program and the finite 
element model of the mast are shown in tables II to IV. In 
these tables the values for the mass moment of inertia and the 
preload vary from zero to their maximum values. The circular 
frequencies from equation (33) and from the finite element 
model of the mast are nearly identical at each mode; thus the 
Fortran program was verified for accuracy. 


7 


Frequency difference, Hz 


Table II. — Comparison of Frequencies 
Obtained From MSC/NASTRAN 
and Equation (33) With 
/ = 0 in. 4 and P = 0 
[M r = 67 lb] 


Modes 

Frequency, Hz 


MSC/NASTRAN 

Equation (33) 

First 

1.27 

1.27 

Second 

9.07 

9.07 

Third 

26.86 

26.87 


Table III. — Comparison of Frequencies 
Obtained From MSC/NASTRAN 
and Equation (33) With 
/ = 0 in. 4 and P = 150 lb 
\M r = 67 lb] 


Modes 

Frequency, Hz 


MSC/NASTRAN 

Equation (33) 

First 

1.25 

1.25 

Second 

8.88 

8.90 

Third 

26.44 

26.71 


Table IV. — Comparison of Frequencies 
Obtained From MSC/NASTRAN 
and Equation (33) With 
/= 38 in. 4 and P = 1501b 
\M r = 67 lb] 



1 2 3 4 5 6 7 8 9 10 11 12 

Mode number 

Figure 7. — Frequency differences from table V as a 
function of mode number. 


Tabic V. — Comparison of Frequencies 
Obtained From Equation (33) for Two 
Cases: / = 0, and / = 38 in. 4 
[P= 150 lb and M { = 67 lb] 



To determine the influence of the mast mass moment of in- 
ertia on circular frequencies, the first 10 circular frequencies 
were extracted and tabulated in table V for 2 mass moments 
of inertia input to the Fortran program, namely, 7 = 0 in. 4 and 
/ = 38 in. 4 The differences between these circular frequencies 
are shown in table V and plotted in figure 7; note that as the 
mode number increases, the influence of the mast mass 
moment of inertia becomes more apparent. The relationship 
between the mode number and the frequency difference in 
figure 7 appears to be that of a second order polynomial. 
MSC/NASTRAN was not included in the study performed 
for table V, since for higher frequencies the mast finite ele- 
ment model demands higher fidelity. This model will need to 
be further refined. 

To incorporate the influence of the mass moment of inertia 
into the finite element model of the mast, the following meth- 
ods were investigated: (1) the mass moment of inertia of 
each BEAM element of the mast was calculated and applied 
to the BEAM property card as the nonstructural mass 
moment of inertia; and (2) the mass moments of inertia of the 
BEAM elements (calculated in step 1) were applied to con- 
centrated mass elements (CONM2), and these elements were 
placed on nodes 1 to 50. As expected, the circular frequencies 
from the finite element model of the mast were not affected at 
all, because the influence of the mass moment of inertia has 
not been included in the Euler beam formulation. 

For a cantilevered beam such as the structure of the mast, 
where most of modal energy is stored in its first few modes, it 
is crucial to predict the first few circular frequencies accu- 
rately. The reason for this is that the dynamic loads are gener- 
ally proportionaj tojhe inverse of the circular frequencies. As 
can be observed in tables II and IV, the frequencies obtained 
from MSC/NASTRAN are accurate. However as the demand 
for higher circular frequencies grows, MSC/NASTRAN 
becomes less accurate. 
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Conclusions 

From the work reported herein, the following conclusions 
may be drawn: 

1. Equation (33) is the eigen condition for the characteris- 
tic equation of the preloaded mast with a tip mass whose 
infinite number of roots are the circular frequencies of the 
structure. 

2. According to our formulations and analyses, MSC/ 
NASTRAN is accurate for predicting the bending circular 
frequencies of masts as long as the influence of the mass 
moment of inertia of the beam is small. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, February 3, 1992 
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